Collective phase synchronization in locally-coupled limit-cycle oscillators 
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We study collective behavior of locally-coupled limit-cycle oscillators with scattered intrinsic fre- 
quencies on d-dimensional lattices. A linear analysis shows that the system should be always desyn- 
chronized up to d = 4. On the other hand, numerical investigation for d = 5 and 6 reveals the 
emergence of the synchronized (ordered) phase via a continuous transition from the fully random 
desynchronized phase. This demonstrates that the lower critical dimension for the phase synchro- 
nization in this system is di — 4. 
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Up to date, much attention has been paid to the collec- 
tive behavior of coupled nonlinear oscillators since those 
systems of oscillators have been known to exhibit re- 
markable phenomena of synchronization Q. Those phe- 
nomena have been observed in a number of physical, bi- 
ological, and chemical systems, and understood rather 
well in terms of various model systems. In the sys- 
tem of globally-coupled oscillators such as the Kuramoto 
model H IK, the mean-field (MF) theory is valid and 
yields analytic results to unveil the phase transition 0. 
Systems of locally-coupled oscillators, on the other hand, 
have not been much studied even though local coupling 
in the system is more realistic in nature. In some existing 
studies 0, 0, E| collective synchronization, in particular, 
frequency entrainment has been investigated. However, 
even numerical results as well as the analytic ones includ- 
ing heuristic arguments do not provide a clear answer 
about the question of the lower critical dimension for the 
frequency entrainment. Phase synchronization has been 
also studied in the previous studies pfl , however, there are 
still many fundamental questions that are not answered. 

In this Letter, we consider a system of locally-coupled 
oscillators on d-dimensional lattices and use the relation 
with a typical model of growing surface, which allows a 
linear analysis to show the absence of synchronization 
up to d = 4. On the other hand, numerical investiga- 
tion performed for d = 5 and 6 reveals the emergence of 
the synchronized (ordered) phase via a continuous tran- 
sition, indicating the lower critical dimension for phase 
synchronization di = 4. 

We begin with the set of equations of motion governing 
the dynamics of N coupled oscillators located at sites of 
a d-dimcnsional hypercubic lattice: 



dt 



(1) 



where fa and u>i stand for the phase and the intrinsic fre- 
quency of the ith oscillator (i = 1, 2, N), respectively. 
The intrinsic frequencies are assumed to be randomly dis- 
tributed according to the Gaussian distribution function 



g{uj) with mean too and variance 2a. For simplicity, we 
set ujq = without loss of generality. The second term 
on the right-hand side represents local interactions be- 
tween the ith oscillator and its nearest neighbors the set 
of which is denoted by Aj . 

Without any interaction (K = 0), each oscillator 
evolves with its own intrinsic frequency, resulting in that 
the system becomes trivially desynchronized. For K > 0, 
the coupling term favors locally ordered (synchronized) 
states and competes against the randomizing force due 
to scattered intrinsic frequencies. When the coupling is 
strong enough to create globally ordered states, the sys- 
tem should exhibit collective synchronization. We here 
focus on phase synchronization which may be probed by 
the conventional phase order parameter 




N 



(2) 



where (• ■ • } denotes the average over realizations of in- 
trinsic frequencies. Phase synchronization is then identi- 
fied by nonzero A in the thermodynamic limit. 

Analytic results are available at the MF level. Namely, 
in the case of globally coupled oscillators, where each 
oscillator is coupled with every other one with equal 
strength K/N, it is known that phase synchronization 
emerges as A ~ (K — K c )@ with (3 = 1/2 near the critical 
coupling strength K c = 2/irg(0) while the correlation 
length diverges as £ ~ \K — K c \~ v with v = 1/2 

When the oscillators are locally coupled, the system 
has been little investigated. Since the nonlinear nature 
of the sine coupling term in Eq. (Qj is the major obsta- 
cle toward analytic treatment, we first suppose that, for 
sufficiently strong coupling strength K, the phase differ- 
ence between any nearest neighboring oscillators is small 
enough to allow the expansion of the sine function in 
the linear regime. With the appropriate continuum limit 
taken in space, the linearized evolution equation for the 
phase </>(x, t) reads 



dt 



w(x) + KV 2 (j) + C(V 4 (; 



(3) 



2 



where oj(x) are uncorrelated random variables, satisfying 
(w(x)) = and (w(x)w(x')) = 2cr<5(x - x') 0. We also 
relax the constraint < <f> < 2ir and extend the range of 
4> to (— oo,oo), for convenience. 

With the irrelevant high order terms neglected, 
this equation reminds us of the celebrated Edwards- 
Wilkinson (EW) equation 10], traditionally describ- 
ing certain surface evolution, by interpreting the phase 
<fi(x, t) as the front height of the growing surface. Note, 
however, that the noise w(x) is generated not by conven- 
tional spatio-temporal disorder but by so-called columnar 
disorder (with spatial dependence only). 

In the context of surface growth models, a central 
quantity of interest is the surface fluctuation width W 
defined by 



w 2 {t) = - d 



[0(x,t) - <£(*)] : 



(4) 



where L is the linear size of the <i-dimensional lattice 
(L d = N) and (j)(t) the spatial average of the phase 
</>(x, t). By means of the Fourier transforms, one can eas- 
ily solve Eq. J2Jl to find in the long time limit (Kt ;§> L 2 
that the steady-state surface width scales for large L 



W 2 



(2a/K 2 )L 4 ' d , d < 4 
(ct/4tt 2 K 2 )]xiL, d = 4 
2(7 /K 2 , d>4. 



(5) 



At any finite values of K, the surface width W thus di- 
verges as L — > oo for d < 4 whereas it remains finite for 
d > 4. This indicates that the surface is always rough 
(except at K — oo) for d < 4 and always smooth (except 
at K = 0) for d > 4. 

It is also straightforward to derive the steady-state 
probability distribution ||: 



P[{0}] ~ exp 



-(K 2 /4a) J (V 2 0)Vx 



(6) 



Notice that the Gaussian property of the probability dis- 
tribution links W analytically to the phase order param- 
eter via A = exp[— W 2 /2}. Therefore our results for W, 
translated into the phase synchronization language, show 
that the oscillators are always desynchronized (A = 0) 
for d < 4 and always synchronized (A ^ 0) for d > 4 in 
this linearized model. 

Our linear theory is valid in the strong-coupling 
regime; as the weak coupling regime is approached, the 
original (nonlinear) system should be more disordered 
than the prediction of the linear theory. This estab- 
lishes that the full nonlinear system described by Eq. 
JTJ should also be desynchronized for d < 4 at any fi- 
nite K. For d > 4, it is reasonable to expect a phase 
synchronization (roughening) transition at a finite value 
of K, although one may not exclude the possibility of 
either the full destruction of the synchronized phase at 



any finite K or the absence of the desynchronized phase 
at any nonzero K. 

Before investigating the full nonlinear system de- 
scribed by Eq. we consider another standard quan- 
tity in surface growth models, the height-height correla- 
tion function C(x, t) = ([</>(x, t) — 0(0, t)] 2 ). In the lin- 
earized regime governed by Eq. 13J), we find the steady- 
state behavior for small x = |x| |8| 



C(x) 



(2<j/K 2 )x 2 L 2 - d , 
(a/2nK 2 )x 2 InL, 
(2a/K 2 )x 4 - d , 



d<2 
d=2 
d>2. 



(7) 



Note that for d < 2 the correlation C(x) diverges with 
system size L, which implies that the average nearest 
neighbor phase (height) difference G = ((V0) 2 } 1 / 2 is un- 
bounded for any finite K in the thermodynamic limit. As 
our linear theory is based on the boundedness of |V</>|, 
there is no range of K where the linear theory applies 
for d < 2. In contrast, for d > 2, G is finite and the 
linear theory is self-consistent at least for large K where 
G(K) < 0(1). We now examine the nonlinear effects due 
to the sine coupling in Eq. Unlike in the linearized 
case, phase <j) may not be bounded even in a finite sys- 
tem but diverge eventually with a finite angular velocity, 
once its intrinsic-frequency term wins over the nearest- 
neighbor coupling term. In the weak-coupling regime (for 
small K), these runaway oscillators with scattered an- 
gular velocities dominate and their phases become com- 
pletely random to one another, leading to the behavior 
A ~ N^ 1 / 2 = L~ d / 2 . On the other hand, in the strong 
coupling regime where the linear theory applies, A van- 
ishes exponentially for d = 3 and algebraically for d = 4, 
with an exponent depending on K [see Eq. 0]. 

We integrate numerically Eq. and measure the 
phase order parameter at various values of K and L for 
d = 2 to 6. For convenience, periodic boundary con- 
ditions have been employed and 2a has been set equal 
to unity. We start from the uniform initial condition 
(4>i = 0) for a given set of {oJi}, chosen randomly accord- 
ing to the Gaussian distribution g(oj) ~ exp(— to 2 /4a), 
and measure the order parameter A averaged over the 
data in the steady state, reached after appropriate tran- 
sient time (Kt ^ L 2 ). Here we have used Heun's 
method to integrate up to 4 x 10 4 time steps, with the 
time step 5t = 0.05, and also average over 100 indepen- 
dent sets of {u>i}. Figure H displays the numerical results 
for the order parameter. For d = 2 and 3, it is clearly ob- 
served that the order parameter decreases rapidly with 
the system size and seemingly approaches zero in the 
thermodynamic limit for any finite K. Detailed finite- 
size analysis shows A ~ L~ d l 2 in the weak coupling 
regime, implying that phases are completely random and 
the system is dominated by runaway oscillators. For 
d = 2, this fully random phase extends to the regime 
of large K, while for d — 3 the linear theory predicting 
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FIG. 1: Behavior of the order parameter A with the coupling 
strength K, plotted in terms of exp(— K), in systems of vari- 
ous size L for (a) d = 2; (b) d — 3; (c) d — 4; and (d) d = 5. 
Symbol sizes correspond to statistical errors of the data. 
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FIG. 2: Log-log plots of A versus L~ for d = 4 at various 
values of K. The data for large K are shown in the inset for 
better visibility. Lines are merely guides to eyes. 



correlated phases [see Eq. Q] appears to work for large 
K, namely, the data fit well to A ~ exp[— [a j ^ K 2 )L\ 
for K > K with G (K ) w 0(1). Numerically, wc find 
that Kq ss ^/2ct/7t 12]. The data for d = 4 seem to 
suggest that for large K A remains finite even in the 
thermodynamic limit, which contradicts our prediction 
based on the linear analysis. To resolve this puzzle, we 
analyze our data carefully by means of finite-size scal- 
ing, and show in Fig. |2Jthe log-log plots of A versus L^ 1 
at various values of K. Manifested for K < 0.28 is the 
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FIG. 3: Log-log plots of A versus L 1 for d 
values of K. Lines are merely guides to eyes. 



5 at various 



fully random phase: A ~ L~ 2 . For K > 0.40, A still 
decreases algebraically with L (see the inset of Fig. EJ : 
A - L~ s(K l It is pleasing that our data for K > 0.40 
agree perfectly with the prediction of the linear theory, 
S(K) = a/8ir 2 K 2 from Eq. (J5J. This result confirms 
that there is no synchronized phase at any finite K for 
d = 4. It would be interesting to explore the possibility 
of a phase transition near K rs Kq = ^/cr/4 between the 
fully random phase and the critical phase described by 
the linear theory; this is currently under investigation. 
For d = 5, it looks evident that there exists an ordered 
(synchronized) phase extended to finite values of K . Sim- 
ilarly to the d = 4 case, the log- log plots of A versus L" 1 
are drawn in Fig. |3| For K < 0.19, we find the fully 
random phase: A ~ L~ 5 / 2 . For K > 0.21, on the other 
hand, A, first decreasing slightly with L, eventually con- 
verges to a non-zero value. In fact, for K > 0.24, this 
saturated value coincides perfectly well with the linear- 
theory value: A = exp[— a/12TT 2 K 2 ). Note here that the 
linear theory breaks down for K < Kq = \J o j% w 0.24 
and the transition into the fully random phase appar- 
ently occurs a little later at K c S3 0.20. It may be very 
interesting to understand this phase transition from the 
stability analysis in the weak coupling limit. 

We next study the critical behavior near the synchro- 
nization transition. In a finite system, we assume the 
finite-size scaling relation 



A = L-P/»f[(K-K c )L 1 /»], 



(8) 



x@ as x 



where the scaling function behaves f(x) 
+oo and f(x) ~ constant as x — > 0. At criticality, it 
leads to 



A(K C ,L) - L-M". 



(9) 



To estimate efficiently the exponent 0/v and the transi- 
tion point K c , wc introduce the effective exponent 

P/v(L) = -HA(L')/A(L)]/HL/L), (10) 

which is expected to approach zero, (3/v, and d/2 for 
K > K c , K = K c , and K < K c , respectively, as L — > oo. 
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FIG. 4: Effective exponent (3/v(L) versus L 1 for d — 5 at 
various values of K. 
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FIG. 5: Data collapse of AL* 3 /" against (_K"/_K" C -1)L 1/ '' in the 
log-log scale for various values of the system size and coupling 
strength. The best collapse is achieved with P/v = 1.4(3) and 
v = 0.45(10). The straight line has the slope 0.63, giving an 
estimation of j3. 

The effective exponent for d = 5, computed at var- 
ious values of K, is plotted in Fig. 21 The data for 
K < 0.19 apparently converge to the weak-coupling value 
5/2, while those for K > 0.21 converge to zero within 
statistical errors. Only the data at K = 0.20 appear to 
converge to a nontrivial value. We thus estimate the crit- 
ical coupling strength K c = 0.200(5) and the exponent 
ratio p/v = 1.6(3). 

To check the finite-size scaling relation directly, we plot 
AI"/" versus (K/K c - l)L l / u in Fig.0and find that the 
data for various values of L and K are best collapsed 
to a curve with choices of K c = 0.200(5), P/v = 1.4(3) 
and v = 0.45(10), which results in p = 0.63(20). As 
expected, the resulting scaling function f(x) converges 
to a constant for small x, and diverges as for large x 
(see Fig. |5J|. 

We summarize our results for d = 5: 

^/^ = 1.5(3), i/ = 0.45(10), K c = 0.200(5). 

Note the apparently substantial deviations from the MF 
values, P/v = 1 and v — 1/2, although the latter may 
not be totally excluded. In view of the argument for the 
MF nature [8j, these apparent deviations are rather un- 



expected and their origin is unclear at this stage. Sim- 
ilarly, we find for d = 6: P/v = 1.0(3), v = 0.45(10), 
K c = 0.158(5), which seem to be consistent with the MF 
values. 

In summary, we have explored the phase synchroniza- 
tion phenomena in the system of locally-coupled oscilla- 
tors with scattered intrinsic frequencies on ci-dimensional 
lattices. A linear analysis shows that the strong cou- 
pling regime can be described by the EW surface growth 
equation with columnar disorder for d > 3. It has been 
shown analytically that the system is always desynchro- 
nized up to d = 4, while numerical integration for d > 5 
has demonstrated the emergence of the synchronized (or- 
dered) phase via a continuous transition from the desyn- 
chronized phase. The lower critical dimension for phase 
synchronization is thus given by di = 4, but the criti- 
cal behavior explored for d = 5 and 6 does not give a 
conclusive result for the upper critical dimension. 
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